ON A MODEL OF MULTIPHASE FLOW 
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Abstract. We consider a hyperbolic system of three conservation laws in one space variable. 
The system is a model for fluid flow allowing phase transitions; in this case the state variables are 
the specific volume, the velocity and the mass density fraction of the vapor in the fluid. For a class 
of initial data having large total variation we prove the global existence of solutions to the Cauchy 
problem. 
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1. Introduction. We consider a model for the one-dimensional flow of an invis- 
cid fluid capable of undergoing phase transitions. Both liquid and vapor phases are 
possible, as well as mixtures of them. In Lagrangian coordinates the model is 
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Here t > and x € K; moreover v > is the specific volume, u the velocity, A the 
mass density fraction of vapor in the fluid. Then A S [0, 1], with A = characterizing 
the liquid and A = 1 the vapor phase; the intermediate values of A model the mix- 
tures of the two pure phases. The pressure is denoted by p = p(v, A); under natural 
assumptions the system is strictly hyperbolic. 

This model is a simplified version of a model proposed by Fan [13] . where also 
viscous and relaxation terms were taken into account. The model is isothermal, see 
(12. ip below; in presence of phase transitions this physical assumption is meaningful for 
retrograde fluids. A study of the Riemann problem for a 2 x 2 relaxation approximation 
of (|l.ip has been done in [10]. We focus here on the global existence of solutions to 
the Cauchy problem for (jl.ip . namely for initial data 

(v,U,X)(0,x) = (v (x),Uo(x),\ (x)) 

having finite total variation. This problem is motivated by the study of more complete 
models, where (jl.l[) is supplemented by source terms. 

The problem of the global existence of solutions to strictly hyperbolic system of 
conservation laws has been studied since long, see [§1 EU [211 US] for general informa- 
tion. If the initial data have small total variation then Glimm theorem [14] applies; 
we refer again to [3] for the analogous results obtained by a wave-front tracking al- 
gorithm as well as for uniqueness and continuous dependence of the solutions on the 
initial data. 

Some special systems allow however initial data with large total variation. For 
the system of isothermal gasdynamics Nishida [19] proved that it is sufficient that 
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the variation TV(v 01 u ) of the initial data is finite in order to have globally defined 
solutions. This result was extended by Nishida and Smoller [20] to any pressure law 
p = k/v 1 , 7 > 1, provided that (7 — l)TV(v ,u ) is small; related results are in [12"] . 
For the full nonisentropic system of 3 x 3 gasdynamics, p = kexp( : ^-s)/v' 1 , for s the 
entropy, Liu [T71 [16] proved the global existence of solutions if (7 — l)TV(v ,p ) is 
small and TV(s Q ) bounded. Temple [26] and Peng [21] obtained similar results. All 
these papers use the Glimm scheme. Analogous results making use of a wave-front 
tracking scheme have been given recently by Asakura [I], [5]; we point out that the 
use of wave-front tracking schemes in case of data with large variation is far from 
being trivial, and a deep analysis of the wave interactions is required. Very general 
results can be proved for systems with coinciding shock and rarefaction curves, [5]; 
however system (jl.ip is not of this type. 

In comparison with the above systems of gasdynamics, in (|1.1[) we keep a 7-law 
for the pressure with 7 = 1, but add a dependence of p on A: we take then p = a(A) jv 
for a suitable function a. System (|1.1[) has close connections to a system introduced 
by Benzoni-Gavage [H] and studied by Peng [22] : it seems however that the proof in 
[2"2"] is not complete. A comparison of these models is done in Subsection 13. 11 We 
mention that also the method of compensated compactness has been applied to (jl.ip , 
see [15], [7] and [HI §12.3, §16], but for different pressure laws. 

In this paper we prove by a wave-front tracking scheme the global existence of 
solutions to (jl.ip for a wide class of initial data with large total variation. We intro- 
duce first a weighted total variation (WTV) of a(A ); this quantity arises in a natural 
way in the problem and has also an analytical meaning, being the logarithmic varia- 
tion in the case of continuous functions. We prescribe a bound on WTV (a(A )); for 
the variation TV(w ,u ) there is not such a bound but, roughly speaking, the larger 
TV(v ,u ) is, the smaller must be WTV (a(A Q )). An important point is that we give 
explicit expressions for these bounds; then our results are qualitatively different from 
some of those quoted above, where a generic smallness is required. 

The plan of the paper is the following. The main result is stated in Section [2l 
Theorem 12.21 The Riemann problem is reviewed in Section [3] together with related 
results; proofs have been given in [2]. The definition of the algorithm is in Section [4] 
The core of the proof are Section [5] - where interactions are studied in detail - and 
Section [5]- where we prove the convergence and consistence of the scheme. A careful 
analysis is needed due to the presence of large waves. 

The paper is completed by two appendices. In the first one we prove the main 
result on the weighted total variation. In the second we study the interaction of 
two shock waves to the light of Section [H namely we look for precise bounds of the 
damping coefficient that controls the reflected wave produced in the interaction; we 
think that this analysis is interesting by its own. Good reading! 

2. Main results. We consider the system of conservation laws (jl.ip . The pres- 
sure is given by 

(2.1) p(v,X) = ^ 

v 

where a is a smooth (C 1 ) function defined on [0, 1] satisfying for every A £ [0, 1] 



(2.2) 



a(A) >0, a'(A)>0, 

2 



see Figure I2TT1 For instance a 2 (A) = ko+X(ki — k ) f° r < fco < k\. As a consequence 
of (|XT]1 and (12~2"1) we have, for every (v, A) G (0, +oo) x [0, 1], 

(2.3) p > 0, p u < 0, p„„ > 0, 

(2.4) p\ > 0, p tlA < . 

Remark that assumptions (|2.3p and (|2.4p are analogous to those usually made on the 
pressure in the full non-isentropic case, [17] . the entropy replacing A. 



< A < 1 




Fig. 2.1: Pressure curves as functions of v. 

We denote U = (v, u, A) G = (0, +oo) x K x [0, 1] and by C/ = (v, u) the projec- 
tion of U onto the plane vu; the same notation applies to curves. Under assumptions 
(|2.1[) and l|2.2p the system (|1.1|) is strictly hyperbolic in the whole SI with eigenvalues 
e i = —\/-Pv(v, A), e 2 = 0, e 3 = A). We write c = j>„ = a(X)/v. 

The eigenvectors associated to the eigenvalues i = 1,2,3, are ri = (l,c, 0), 
r 2 = (—p\,0,p v ), r 3 = (— l,c, 0). Because of the third inequality in (|2.3|) the eigen- 
values ei, e 3 are genuinely nonlinear with Ve^ • r^ = p vv /(2c) > 0, z = 1,3, while 
62 is linearly degenerate. Pairs of Riemann invariants are R\ = {u — a(A) log v, A}, 
R 2 = {u,p}, R 3 = {u + a(A) logv, A}. 

We denote by TV(/) the total variation of a function /. In the case / : R — > 
(0, +oo) we define the weighted total variation of / by 



WTV(/) = 2sup^ 



1/(^0-/(^-1)1 



where the supremum is taken over all n > 1 and (n + l)-tuples of points Xj with 
x < xi < . . . < x n . This variation is motivated by the definition (|3.6p of strength 
for the waves of the second family. If / is bounded and bounded away from zero then 
clearly 

' TV(/) < WTV(/) < r^TVf/) . 



sup/ inf/ 
PROPOSITION 2.1. Consider / : K — > (0, +oo); i/ien 

(2.5) ^4 TV (log(/)) < WTV(/) < TV (log(/)) . 

sup/ 



Moreover, if f £ C(R) tfien WTV(/) = TV(log(/)). 
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The proof is deferred to Appendix [X] In (|2.5[) . in the inequality on the right, the 
strict sign may occur if / is discontinuous; see Remark I A. II 

We provide system with initial data 

(2.6) U(x,0) = U (x) = (v (x),u (x),X (x)) 

for i£l Denote a (x) = a (X (x)), p (%) = p (v (x), X (x)); remark that inf a {x) > 
a(0) > 0. The main result of this paper now follows. 

Theorem 2.2. Assume 112.1]) , 112.2]) . Consider initial data 112.6]) with v (x) > 
v > for some constant v and < X {x) < 1. For every m > and a suitable 
function k(m) £ (0, 1/2) the following holds. If 

(2.7) TV (log(p„)) + — !— TV(u ) < 2(l - 2WTV(a„))m 

mt a Q V / 

(2.8) WTV(a Q ) < k{m) 

then the Cauchy problem Sl.l]) , \2.0jl has a weak entropic solution (v,u,X) defined for 
t € [0,+oo). Moreover the solution is valued in a compact set of and there is a 
constant C(m) such that for every t € [0, +00) 

(2.9) TV(v(t,-),u(t,-))<C(m). 



The function k(m), whose expression is given in (|6.2ip . deserves some comments. 
The interaction of two waves a, a' of the same family i = 1, 3 produces a wave (3 of 
the same family i and a "reflected" wave 5 of the other family j (J = 1, 3, j ^ i). For 
a suitable definition of the strengths of the waves we prove that \S\ < d-mm{\a\, \ at'\} 
for a damping coefficient d < 1 depending on a and a', see Lemma 15.61 The function 
k above depends essentially on the supremum of such coefficients d; we prove that 
fc(0) = 1/2 and that k(m) decreases to as m — > +00. In particular then WTV(a G ) < 
1/2. The assumptions (|2.7|) . (|2.8|) read as analogous to those in [20]: the larger is m, 
the smaller is k(m), and vice- versa. The occurring of a possible blow-up when the 
bound on WTV(a Q ) does not hold is an interesting open problem. 

The variation of A G appears both in condition (|2.7p . because of p Dl and in (|2.8p . 
Using the definition of the pressure, we can replace (12. 7|) by the slightly stronger 
condition TVlog(w ) + 2TVlog(a ) + -^TV(u„) < 2(1 - 2WTV(a Q ))m or even 

1 / 2m +1 \ 

TVlog( Wo ) + ^TV( Uo ) < 2m 1 TV(log(a )) 

a((J) \ m J 

by making use of (|2.5p . In particular if A is constant we recover the famous result 
by Nishida [19]. 

Clearly X(t,x) — X (x) for any t because of the third equation in (jl.ip ; this is 
why only v and u appear in the estimate (|2.9p . In other words system (II. ip can be 
rewritten as a p-system of two conservation laws with flux depending on x, namely 
for the pressure law p — p(v, X (x)) — a 2 (A (a;)) /v. 

The proof of Theorem 12.21 makes use of a wave-front tracking scheme where we 
exploit the special structure of system (II. ip by differentiating the treatment of 1 and 
3 waves from that of 2 waves. Our algorithm is a natural extension of that in [3], 
where the system for A G constant is studied, in presence of a relaxation term. 
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Here we consider a linear functional as in [3] that accounts for the strengths of 
all 1 and 3 waves, with a weight £ > 1 assigned to shock waves; a crucial point in 
the proof is the choice of £ as a function of m. This functional differs with that in 
[19] . [4], where £ is missing and only the variation of shocks is taken into account. 
Moreover, motivated again by [3], we do not introduce a simplified Riemann solver for 
interactions between 1 and 3 waves but only for interactions involving the 2-contact 
discontinuities. The interaction potential considers then uniquely interactions of 2 
waves with 1 or 3 waves approaching to it. 

System can be written in in Eulerian coordinates. Denoting p = 1/v the 

density, the pressure law becomes p — a (X)p and (jl.ip turns into 

!Pt + (pu) x = 

(pu) t + (pu 2 + p(p, X)) x =0 
{p\)t + {p\u) x = . 

A global existence result of weak solutions for (|2.10|) holds by Theorem 12.21 because 
of [271. 

3. Preliminaries. 

3.1. Comparison with other models. In 6J many models for diphasic flows 
are proposed and studied. In a simple case (no source terms, the fluid in either a 
dispersed or separated configuration) and keeping notations as in [6J page 35], they 
can be written as 

r {piRi) t + {piRiui) x =0 

(3.1) I {p g R g ) t + {p g R g u g ) x =0 

[ {piRlUl+ PgRgUg)t + {piRluf + Pg RgU 2 g +p) X =0. 

Here the indexes I and g stand for liquid and gas. Therefore pi, Ri, ui are the liquid 
density, phase fraction, velocity, and analogously for the gas; clearly Ri + R g = 1. The 
pressure law is p — a 2 p g , for a > a constant. Equations (|3.1[) state the conservation 
of mass of either phases and the total momentum. 

A case studied in [6l page 44] is when ui — u g and pi is constant, say equal to 1. 
The unknown variables are then i?/, u, p g , and it is assumed < Ri < 1 and p g > 0; 
as a consequence < R g < 1. Under these conditions, and writing still pi instead of 
1 for clarity, we define the concentration c = > and deduce the pressure law 

P = Q 2 c 1 ^j ii . We obtain exactly the model of [2"2"] : 

( {Ri) t + {Riu) x =0 

(3.2) <^ (R l c) t + {R l cu) x =0 
{ {R l (l + c)u) t + (Ri(l + c)u 2 +p) x =0. 

This system is strictly hyperbolic for c > 0. Remark that the three eigenvalues of 
(13. 2p coincide with u at c = and if c vanishes identically then (|3.2[) reduces to the 
pressureless gasdynamics system. System (|3.2[) is analogous to (12.10[) but the pressure 
laws are different. In fact the variables p and A of (|2.10| write p — piRi + p g R g and 

A = Pl Ri+p g R g = TTc and then R l = i l ~ Pg = p-i-\l-X) - If we sum U P the tw0 

first equations in (|3.2p we find the first equation in (I2.10p : the third (resp. second) 
equation in (|3.2p becomes the second (third) equation in p.lOp . The choice p = a 2 p g 
for the pressure in (|3.ip gives p = a 2 i_^_ x ^ ■ 
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Notice that the pressure vanishes in presence of a pure liquid phase and this is 
the main difference with (|2.ip . (|2.2[) . 

We compare now (|2.10[) and (|3.2|) in Lagrangian coordinates. Consider for (|3.2p 
the change of coordinates y — Ridx — Riudt based on the streamlines of the liquid 
particles (because Ri = piRi), [22]. Denote w = — 1 = ^ = — . Then for p = ^ 



Pa 

system (13.21) turns into 

Wt — u x =0 

(3.3) { ({l + c)u) t + Px =0 

c t =0. 

It is more interesting however to consider for system (|3.2I) the change y — (l+c)Ridx— 
(1 + c)Riudt = pdx — pudt into Lagrangian coordinates based on the streamlines of 
the full density p. Let w be as above and v = = -y-. Then system (|3.2j) 
becomes system (jl.ip with a 2 (A) = a 2 jf^ = a 2 A. As a consequence the pressure 
law p(v,X) ~ a 2 (X)/v does not satisfy (|2.2p . This difficulty can be overcome as 

2 

follows. Fix any < a\ < a,2 < a and consider for Cj = -tzs the invariant domain 
{(Ri,u, c): < ci < c < c 2 }, [12]. In this domain < bi < A < 6 2 < 1, for h = a 2 /a 2 . 
If we denote /i = then the function 6(^i) = a(A) = a(6i + (62 — &i)/i) makes the 

pressure law p(v, A) = b 2 (p)/v, with G [0, 1], satisfy both conditions in ()2.2|1 . 

3.2. Wave curves and the Riemann problem. In this section we recall some 
results about the wave curves for system and the solution to the Riemann prob- 
lem; see [2] for more details. 

The shock-rarefaction curves through the point U = (v 0l u , X ) for are 

(3.4) $i{v,U ) = (v,cj )l (v,U ),X ) , » = 1,3 

u a + a(X ) ■ (v — v )/ y/vv Q v < v shock 
u Q + a(X ) log(v/v ) v > v , rarefaction, 



4>i (v, U ) 
&3(v, U ) 



u Q — a(A G ) log(v/v ) v < v rarefaction 

u a — a(X ) ■ (v — v )/ y/vv v > v , shock, 



(3.5) < &2(A, U ) — (v \ ,u ,X) , A e [0, 1] contact discontinuity. 
V a ( A °) / 

The curves $1, $2 and $3 are pZane curves: $1 and $3 lie on the plane A = A Q 
while $2 on u = u a . 

Definition 3.1 (Wave strengths). Under the notations {3.$ , ( TO)) we define the 
strength Si of a i-wave as 

1, f v\ a(X)-a(X ) 1 (v 

3.6 £l = -log — , e 2 = 2 £ 3 = logf — 

2 V w o/ a(A) + a(A ) 2 \v 

According to this definition, rarefaction waves have positive strengths and shock waves 
have negative strengths. Given the initial datum A Q = X Q (x), denote 



(3.7) a* = supa(A (x)) , a* = inf a(A Q (a;)) 



r 1 • a 

a * = — 
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e if £ > , 

sinh e if e < . 



Then [a]* < a(ij+a(o) < ^ anc ^ l £2 l — 2 [a]* < 2. It is useful to define also the function, 
see [22], 

(3.8) h(e) -- 
Then we have for i = 1 , 3 

(3.9) <t>i(v, U ) = u + a(X ) ■ 2h(si) . 

At last we consider the Riemann problem. This is the initial-value problem for 
(jl.ip under the piecewise constant initial condition 



(3.10) (v,u,X)(0,x) 



(vt, ut, X t ) = U e if x < 
(v r ,u r ,X r ) = U r if x > 



for Ui and U r in CI. We denote a r = a(A r ), p r = a^,/v r , and similarly a^, p^. 

Proposition 3.2. Fix any pair of states Ui, U r in Cl; then the Riemann problem 
\3. 1 0\) has a unique Cl-valued solution in the class of solutions consisting of 
simple Lax waves. If Si is the strength of the i-wave, i — 1,2,3, then 

£ 3 - £i = ^ log i 2{ath{e\) + a r h(e 3 )) = u r - u e . 

Moreover, let v > 6e a fixed number. There exists a constant C\ > depending on 
v and a(A) swc/i that if Ui, U r € CI and v > v, then 

(3.11) leil + M + N^CilE^-E/rl. 



For the proof, see [2]. One can easily find that 

(3.12) N + M < i|log(p r )-lo g ^)| + -— — I T \u r -u e \ 

< \\ logK) - log(^)| + | log(a r ) - log(a £ )| + . ] f \u r - u £ \ . 

2 2 mm{ ai,a r \ 

We remark that for any Riemann data (vi,U£, Xi), (v r ,u r ,X r ), the A component of 
the solution takes value A^ for x < and A r for x > 0. The fact that the interfaces 
between different phases are connected by a stationary wave can be interpreted then 
as a "kinetic condition" , pQ , analogous to Maxwell's rule. 

4. The approximate solution. In this section we define a wave-front track- 
ing scheme [5] to build up piecewise constant approximate solutions to p.l[) . More 
precisely we follow the algorithm introduced in [3]. 

First, we approximate the initial data. For any v € N we take a sequence 
{v v Q , Uq, Aq) of piecewise constant functions with a finite number of jumps such that 

(i) TVp" < TV Poi TV< < TVu , WTVa(A^) < WTVa(A D ), inf a" > inf a ; 

(ii) Um a: ^._ 0O (^,^,A^)(a;) =Ym\ x ^- 00 {v ,u ,X )(x); 

(iii) ||(<X,Ao)-K^o,A )|| L x<± 
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where p v = a 2 (X^)/v^. Second, we define the approximate Riemann solver. We 
introduce positive parameters r\ = r\ v , p = p v \ they control respectively the size of 
rarefactions and the threshold when a simplified Riemann solver is used. Define also a 
parameter s > strictly larger than all possible speeds of wave-fronts of both families 
1 and 3. These parameters will be determined at the end of Section [6] 

• At time t = we solve the Riemann problems at each point of jump of 
(Vg, Ug, Aq)(0+, ■ ) as follows: shocks are not modified while rarefactions are 
approximated by fans of waves, each of them having size less than r\. More 
precisely, a rarefaction of size e is approximated by N = [e/77] + 1 waves whose 
size is e/N < 77; we set their speeds to be equal to the characteristic speed of 
the state at the right. 

Then (v, u, A) (t, •) is defined until some wave fronts interact; by slightly chang- 
ing the speed of some waves [9] we can assume that only two fronts interact 
at a time. 

• When two wave fronts of families either 1 or 3 interact we solve the Riemann 
problem at the interaction point. If one of the incoming waves is a rarefaction, 
after the interaction it is prolonged (if it still exists) as a single discontinuity 
with speed equal to the characteristic speed of the state at the right. If a new 
rarefaction is generated, we employ the Riemann solver described before and 
divide it into a fan of waves having size less than 77. 

• When a wave front either of family either 1 or 3 interacts with a 2-wave we 
proceed as follows. Let 82 be the size of the 2-wave and 6 the size of the other 
wave. 

— If \5i6\ > p we solve the Riemann problem as above, that is with the 
accurate Riemann solver. 

— If |$2<J| < p we prolong the 1- or 3- wave with a wave of the same family 
and size. Since the two waves do not commute, a non-physical front 
is introduced, [5], with fixed speed s > 0. The size of a non-physical 
wave is set to be |u r — U(\, where ui, u r are the u components of the left 
and right states of the wave. We call this solver the simplified Riemann 
solver. 

• When a non-physical front interacts with a front of family 1, 2 or 3 ("phys- 
ical"), we prolong the solution with a physical wave of the same size and a 
non-physical one, computing the intermediate value consequently. 

We refer for the last two items to Proposition 15. 121 below. Remark that two non- 
physical front cannot interact since they have the same constant speed s. We denote 
by MV the set of non-physical waves. 

5. Interactions. Fix the index v introduced in the previous section. We shall 
prove in Subsection 16.11 that the algorithm described above is defined for any t > 
and provides for any initial data (v", u v a , A£) a piecewise constant approximate solution 
(v u , u 11 , A") = (v,u, A), where we dropped for simplicity the index v. Here we study 
the interaction of waves. 

For K np > and t > we define the functional L and the interaction potential 
Q, both referred to (v,u, X)(t, •), by 

L{t)= 2J + K np L np , L np = ^2 M 

i=l,3 j&MP 

(5.1) Q(t)= 1*1+ E MM- 

73 at the left of 82 71 at the right of 82 
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Remark that L takes only into account the strengths of both 1 and 3 waves and that 
of non-physical waves. For contact discontinuities we define 

i cd = ^|7 2 |=WTVa(A^). 

Finally, for £ > 1 and K > we introduce 

(5.2) Lc — -^rarefactions £ -^shocks H~ K n pL n p , 

(5.3) F = L^ + KQ . 

For simplicity we omitted to note the dependence on K np in the functional and on 
K np , £, K in F; the choice of K np shall depend on that of K, see Proposition 15. 121 

Observe that, if A is constant, then Q — and F = L^, whose variation was 
analyzed in Lemma 3.2 of [3]. Hence we will assume from now on that 

(5.4) A = WTV(a Q ) > . 

By assumption (i) in Section [4J one has L cc i < A Q . 

In the following sections we analyze in detail the different types of interactions. 
Recalling the definition of h, (|3.8|) . and with the notation of Figure IS~T1 we introduce 
the following identities, see (3.1), (3.2) in [22] : 

(5.5) e 3 - ei = a 3 + fa - ai - j3i 

(5.6) a e h(ei) + a r h(e 3 ) = ath{a\) + a m h(a 3 ) + a m h(/3i) + a r h(j3 3 ) . 




Fig. 5.1: A general interaction pattern. 

Formula (|5.5[) does not depend on A and follows easily by equating the specific volumes 
v before and after the interaction time. By equating the velocities u we obtain (|5.6p . 
These properties are a consequence of the definition (|3 - 6[) of the strengths for 1- and 
3- waves and of (|3.9[) . 

5.1. Interactions with a 2-wave. We consider first the interactions of 1 or 3 
waves with a 2 wave, see Figure I5T21 

PROPOSITION 5.1 (0). Denote by X e , A r the side states of a 2-wave. The 
interactions of 1- or 3-waves with the 2-wave give rise to the following pattern of 
solutions: 



interaction 


outcome 




Xi 


< A r 


Xi > A r 


2 x IR 


1R + 


2 + 3R 


li? + 2 + 35 


2 x 15 


15 + 


2 + 35 


15 + 2 + 3i? 


3R x 2 


15 + 


2 + 3i? 


IR + 2 + 3R 


35 x 2 


LR + 


2 + 35 


15 + 2 + 35. 
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The next lemma is concerned instead with the strengths of waves involved in the 
interaction above. The inequalities (|5.8|) improve the inequality (3.3) in [22) in the 
special case of two interacting wave fronts, one of them being of the second family. 
More precisely under the notations of L 22j we find a term l/(a r + at) instead of 
1/ min{a r , ai}. The proof differs from Peng's. Our estimates are sharp: in some cases 
(|5.8jl reduces to an identity. 





Fig. 5.2: Interactions, (a): from the right; (6): from the left. 

Lemma 5.2 ([2J). Assume that a 1 wave of strength 5i or a 3 wave of strength 
63 interacts with a 2 wave of strength 82 — 2(a r — ae)/(a r + at). Then, the strengths 
Si of the outgoing waves satisfy: £2 = 82 ctnd 



(5.7) 



\e i -S i \ = \£ j \<-\6 2 \-\6 i \ < [a] m \5i\ 



for i,j = 1,3, % 7^ j. Moreover, 



(5.8) 



£1 



\S S \ + \S 3 \[S 2 ] 



if 1 interacts 
if 3 interacts . 



Here [x] + = max{x,0}, [x]_ = max{— x, 0}, x € M. Remark that the colliding 
1 or 3 wave does not change sign across the interaction. Moreover the functional L 
increases if and only if the incoming and the reflected waves are of the same type; this 
happens when the colliding wave is moving toward a more liquid phase. 

Now we prove that F is decreasing for suitable K when an interaction with a 
2-wave occurs. The potential Q is needed to balance the possible increase of Lj. 

Proposition 5.3. Assume A a < 2 and consider an interaction of a I or 3 wave 
with a 2 wave, with the notation of Lemma \5.2[ Then AQ < 0. If moreover 



(5.9) 
then 
(5.10) 



K > 



2£ 



A n 



<f |AQ| 



and hence AF < 0. 

Proof. We consider the interaction of a 3-wave with a 2-wave, as in the proof 
of Lemma 15.21 see [5] ; the symmetric case follows in an analogous way. We use the 
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notation as in Figure IST27 a). We define L* d = L~ d + L+ d , L^ d meaning right or left of 
the 2 wave under consideration. 
By assumption, one has 



(5-11) L cd = L- d + L+ +\S 2 \= L* cd + \5 2 \ < A < 2 . 

Recall that e\ — 5% = £3 and 

(5.12) | £l |-|<5i| = M, if 5 2 > , 

(5.13) |^| = -leg], if 5 2 < , 



so that in particular | ^3 1 = — |<5i||. An estimate for AQ follows at once because 
of (ETTT1) : 

AQ=-|^ 1 | + (| ei |-|«J 1 |)L- d + | £3 |£+ d < \\5M{Ll d -2) 
(5.14) < i|d 2 * 1 |(A,-2) <0. 

Hence, using (|5.7[) . we get 

(5-15) ekal + f AQ < i|^ x | L + - 2)1 < 

because of (|5.9p : this proves (|5.10p . Finally, by using (|5.15[) we get 
(5.16) AF = AL^ + KAQ < £|e 3 | + - 8 X \ + KAQ < 0. 

□ 

5.2. Interactions between 1 and 3 waves. Here we analyze the possible in- 
teractions between 1- and 3- waves. Two situations may occur, see Figure HT3l either 
the waves belong to different families or they both belong to the same family. In this 
last case, at least one of the waves must be a shock. 



£1 £3 £1 £3 




Fig. 5.3: Interactions of 1 and 3 waves. 

Lemma 5.4 (Different families interacting). If a wave of the third family interacts 
with a wave of the first family, they cross each other without changing their strength. 

Proof. See also Lemma 3.1 in 3 . Using notation as in Figure l5~37 a) we have 
£3 — £1 = ^3 — 8\ and h{e\) + ^(£3) = h(Si) + h(Sz). The uniqueness of solutions to 
the Riemann problem implies E\ = Si, £3 — 83. 

Remark that here AL^ = = AQ and then AF = for all £ > 1 and K. D 

Lemma 5.5 (Same family interacting: outcome). Assume that a wave 03 of the 
third family interacts with a wave (3% of the third family, giving rise to waves E\, £3. 
Then 
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(i) a 3 < 0, (3 3 < =>■ £i > 0, e 3 < 0, 

(ii) a 3 f3 3 < =*> ei < 0. 

An analogous result holds for interacting waves of the first family. 

Proof. The proof can be done in a geometric way by observing the mutual posi- 
tions of the curves, 19, 25J. A simple alternative proof by analytical arguments now 
follows. We have 

(5.17) £ 3 - £i = Q 3 + (3 3 

(5.18) fc(ei) + h(e 3 ) = h(a 3 ) + h{[3 3 ) . 

In case (i) these formulas read E\ — £ 3 = |a 3 | + |/3 3 | > and —h{e\) — h(e 3 ) = 
sinh(|a 3 |) + sinh(|/? 3 |) > 0. If it were £ 3 > then e\ > from the first equality and 
£i < from the second, a contradiction. Therefore £ 3 < so that £i + |e 3 | = |a 3 | + \(3 3 \ 
and — h(ei) + sinh(|e 3 |) = sinh(|a 3 |) +sinh(|/3 3 |). Analogously, if it were e\ < 0, using 
elementary inequalities we get = sinh(|ei|) + sinh(|a 3 | + \(3 3 \ + |ei|) — sinh(|a 3 |) — 
sinh(|/3 3 |) > 2sinh(|ei|), a contradiction again. Hence e\ > 0. 

In case (ii) assume a 3 < 0, (3 3 > 0; the other case is dealt analogously since (|5.17p . 
(|5.18|) are symmetric in a 3 , /3 3 . We have e 3 — e± — —\ota\ + \P 3 \ and h(e\) + h{e 3 ) — 
-sinh(|a 3 |) + |/3 3 |, then [h(ei) + e 1 } + [h(e 3 )~e 3 } = \a 3 \ -sinh(|a 3 |) < 0. Ife 3 > 0, this 
last equality becomes h{e\) +£\ = \a. 3 \ — sinh(|a 3 |) < that implies E\ < 0. If e 3 < 0, 
then h(e\) +£\ = [\a 3 \ — sinh(|a 3 |)] — [|e 3 | — sinh(|e 3 |)]. If it were S\ > it would be 
l a 3| < | £3 1) since the map x i— > x — smhx is decreasing; but from \e 3 \ + \e\\ = \a 3 \ — |/? 3 | 
we would get that |e 3 | < |a 3 |, a contradiction. Hence in all cases one has e\ < 0. □ 

Now we give sharper estimates for the interaction of waves of the same family: we 
prove that the strength of the reflected wave is bounded by the size of each incoming 
wave, multiplied by a damping factor smaller than 1 . This property will be crucial in 
the next section and it holds also for interactions with a 2-wave, with damping factor 
[a]*, see (|5.7I) . In the case below however the coefficient depends on the strengths of 
the incoming waves; this happens also when non-physical waves are generated, see 
Proposition 15. 121 We assume that 

„s the strength of any interacting i-wave is less than m, 

for some m > and i = 1,3. 

In the special case of interaction of waves of the same family producing two outgoing 
shocks we give a more precise result in Appendix [B] 

Lemma 5.6 (Same family interacting). Consider the interaction of two waves of 
the same family, of sizes on and Pi, i — 1,3, producing two outgoing waves £\, e 3 ; 
assume H5.19\) . Then the following holds. 

(i) There exists a damping coefficient d — d(m), with < d < 1, such that 

(5.20) \sj\ <d(m)-min{|oi|,|A|}, j^i. 

(ii) // the incoming waves are both shocks, the resulting shock satisfies \si\ > 
max{|ai|, If the incoming waves have different signs, both the amount of shocks 
and the amount of rarefactions of the i th family decrease across the interaction. 

In any case 

(5.21) |e»| < k| + |&|. 
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Proof, (i) To fix the ideas, assume i = 3. We have 

(5.22) e 3 - e x = a 3 + ft 

(5.23) h( El ) + h{e 3 ) = h(a 3 ) + h({3 3 ) 

and then 

(5.24) h( Sl ) + h{ex + a 3 + ft) = h{a 3 ) + h{(3 3 ) . 

Remark that the equation (|5.24|) is symmetric in a 3 , ft; from the implicit function 
theorem we find that E\ = £i(a 3 ,ft) is C 1 . Using the notation E\ = r, 0:3 = a, 
ft = b, the identity (|5.24|) rewrites as 

(5.25) h(r) + h( T + a + b)- h(a) - h{b) = , 
with t = r(a, b). One verifies that r(a, 0) = r(0, 6) = and that 

_ h'(a)-h'(T + a + b) _ h'(b) - ti(T + a + b) 

Ta ~ h'(r) + h'{r + a + b)' n ~ h'(r) +h / {r + a + b)' 

As a — > 0, one has r a — >• (1 — h'(b))/(l + h'(b))\ then |t o (0,6)| < 1 and it can be 
bounded by a positive constant less than 1 that depends on m. The same argument 
works for 77,. 

To complete the proof, we show that \t\ < min{|et|, in the non-trivial case 
a 7^ b. We argue by contradiction, using an argument of [53]. Suppose that 
|t| > \a\; we can assume r > 0, since the case r < can be proved by using the 
equality (|5.25p written in terms of Git) — —h(—t). Since the function h is increasing 
we have, for r > |o|, 

h(r)>h(a), h(r + a + b) > h(b) 

Moreover, one of the two inequalities is strict: if a < the first, if a > the second. 
Hence we contradict (|5.25p . 

(ii) From [3J we already know that AL = \e\\ + \e 3 \ — \a 3 \ — \(3 3 \ < , hence (|5.21[) . 
If the incoming waves are both shocks, then (|5.22p becomes 

(5.26) |£3| + |ei| = |a 3 | + |/3 3 |. 

From (i) we have |ei| < \a 3 \, \(3 3 \ and hence the first part of (ii). On the other hand, 
if a 3 p 3 < 0, we have E\ < 0. We can assume a 3 < < (3 3 ; hence (|5.22[) becomes 
£3 = I /3 3 1 — |« 3 1 — \e\\. If £3 > 0, then |e 3 | < |/3 3 |; if £3 < 0, using (i) again one finds 
M < |a 3 |. □ 

Remark 5.7. The damping coefficient d(m), see Figure I5T41 is given by 

d(m) = max . ,. . ,, n , 
l;i|»min{|a|,|6|} 

where the function s(a, b) satisfies h(e) + h(e + a + b) — h(a) — h(b) — 0, see (|5.24[) . 
Hence d(m) increases with m, and vanishes as m — > because quadratic interaction 
estimates hold for m small. 

Moreover, it is asymptotic to 1 for m large. Indeed, from the proof of Lemma [STB] 
we have r a (0,b) = i^g; then r a (0,b) = if b > and |r a (0,6)| = ^Mg_i < 
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Fig. 5.4: The coefficient d(m) 



c°shm+i if b < °- Therefore |r a (0,6)| < coshm+i for ever y & > and an analogous 

cosh m— 1 
cosh rn+1 



estimate holds for \r b (0,a)\. Hence d(m) > £2 ^f^r = c(m); we refer to Lemma TB. II 



for the role of this quantity. 

Remark 5.8. When a rarefaction interacts with a 1- or 3-wave, its size does not 
increase. Indeed, the size does not change upon interactions with waves of the other 
family by Lemma l5.41 if the rarefaction interacts with a shock of the same family we 
apply Lemma \5.6V ii). Remark moreover that by Lemmas 15.41 [531 a rarefaction never 
produces a rarefaction of the other family by interactions with 1- and 3- waves. 

Remark 5.9. If two waves of the same family interact, the wave belonging to 
that family can be missing, while the "reflected" wave is always present. This follows 
easily from (|5~25|) , 1(525)1 . 

Proposition 5.10 (Variation of F). Consider the interactions of any two wave 
fronts of the same family, 1 or 3, and assume H5.19\) . If 

(5.27) 1 < £ < - and K < ^ ~ ' 



d A 
then AL C < and AF < . 

Proof. Let two waves a,, /Jj interact, i = 1,3, giving rise to waves £i, £3. We 
consider i — 3, the other case being analogous. Using (|5.2ip . we get 

(5.28) AQ = (| E3 |-|a 3 |-| J 9 9 |)L++| ei |L^<| ei |Ii^ < MA,. 
Now we claim that 

(5.29) AL C + |ei|(£- 1) < 0. 

From this estimate it follows AF = AL^ + KAQ < \ei\ (1 - f + KA Q ) < because 
of (|5.27p . To prove our claim we consider the possible cases; we make use of (|5.22l) . 



SS -> RS Since AL = then 



(5.30) + (^ - l)|ei| = + N - |a 3 | - |/3 3 |) < 0. 



SR, RS 5ii Assume a 3 < < #3; then |5T29l reads (2£ - l)|ei| + |e 3 | 



1/^3 1 — ^l a 3| ^ 0- For later use we prove the stronger inequality 
(5.31) e\£i\ + \e 3 \-\f3 3 \-£\a 3 \ <0. 
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Indeed, from Lemma we have |e 3 | < |/3 3 |, while £|ei| < |as| from (|5.20|) . (|5.27[h . 
Assume a 3 < < [3 3 ; then (|5T29|) is (2f-l)|ei| 



SR, RS — ► SS 



-£(M 



"3 



\0a\ < 0. We prove also in this case the stronger inequality 
(5.32) f a| ei | +c( | e3 |_| Q3 | ) _| /% |< _ 

Indeed, by (|5.17p and again because of (|5.20|) , (|5.27[) i , one has 

£ 2 N + £(M - \a 3 \) - \0 3 \ = e\si\ + CCkil - Iftl) - Iftl 

= (C + l)(C| ei |-|/%|)<0. 
This proves the claim and concludes the proof. □ 

Remark 5.11. From the above proof we see that AL% < for £ = 1. This was 
a key point in [T!5], where however a different choice of strengths was done. In [3] the 
inequality AL^ < was proved to hold also for 1 < £ < £ Q , for some £ Q > 1; the 
condition (|5.27P i gives an estimate of such a threshold. 

More precisely, in the first two cases of Proposition 15.101 we have AL^ < for 
every £ > 1. The third case is analyzed in detail in Lemma TB. II we prove there that 
AL 5 < for any £ > 1 if c(m) < 1/2, while we need 1 < £ < 2c (m)-i ^ C ( TO ) > V 2 - 

5.3. Non-physical waves. In this subsection we compute the strength of a 
non-physical wave generated by an interaction and prove that it does not change in 
subsequent interactions. We introduce the following notation: given Ui = (vi,ut,Xi) 
and A r we define by 

U* lr = $ 2 (A r , U t ) = {A r£ v e , u e ,X r ) 



the state on the right of a 2-wave with left state Ut 
right, where A r g = a 2 {X r ) / a 2 {Xg) . See (|3~5]) and [2]. 



(v£,U£,Xi) and A = X r on the 





Fig. 5.5: Simplified Riemann solver 

Proposition 5.12 (Non-physical waves). Consider Ui — (v£,ug,\i). Let U r = 
(v r , u r , X r ) be connected to U£ r by a 1-wave of size Si and U q = (v q , u q , Xe) be connected 
to Ue by a 1-wave of size 5±, see Fiaure \575\( a) . Assume 15.19\) . 

Then U* r and U r differ only in the u component; if S2 denotes the size of the 
2-wave, there exists a constant C Q — C (m) such that 



(5.33) 



\\U* r -U r \\ = \u q -u r \ <0,Mi| 



A similar result holds for the interaction of a 3-wave, see Fiaure ld^oY b ) , again under 

fosp . 
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Moreover the size of a non-physical wave does not change in subsequent interac- 
tions. For any K > and K np < K/C a at any interaction involving a non-physical 
wave we have 

(5.34) AF<0, 

with AF < when a non-physical wave is generated. 

Proof. Recalling [2j Lemma 2], only the u component will be different after 
commutation of the 1- and the 2-wave. We find that 



hence 



u q —ui = 2aih{5\), u r — ug = 2a r h(5\) 

smh rj 

\u q — u r \ = 2\a£ — a r \ ■ \h(5i)\ < \S 2 Si\ ■ 2a(l) max . 

(;•• n- m rj 

Then ([5331 follows with C (m) = 2o(l) • 

Next, assume that a non-physical wave interacts with a 2-wave. Since the values 
of u do not change across a 2-wave, the left and right values of u of the non-physical 
wave do not change across the interaction; hence the size does not change. 

Assume then that a non-physical wave interacts with a 1- or 3- wave of size 5. Since 
A is constant, we refer only to the components v, u. Let (ye,U() and (vg,u q ) be the 
side states of the non-physical wave before the interaction and (vt, u q ), (v r , u r ) be the 
side states of the physical wave. After the interaction, let (vi,ui) be the intermediate 
state. One has 

u r — u q = 2a(X)h(5) — ug — U£ , 

hence \u q — ui\ = \u r — ut\. 

At last, we consider the functional F. The potential Q is unaltered when non- 
physical waves interact with other waves. The only cases in which L np changes are 
when a non-physical wave arises. Assume that a 1- or a 3-wave of size 8 interacts with 
a 2-wave of size 82, producing a wave of same size and a non-physical wave. Then 
AQ = -\S8 2 \ and AL 5 = K np AL np < K np C \S5 2 \; hence AF = AL e + KAQ < 
\5S 2 \(K np C - K), then (15341) . □ 

5.4. Decreasing of the functional F and control of the variations. We 

collect first the previous results into a single proposition. 

Proposition 5.13 (Local decreasing). Consider the interaction of any two waves 
either of families 1, 2, 3 or non-physical. Assume i f 5. 1 9\) for some m > 0; let C a — 
C {m) as in Proposition \5. 1£[ Finally let A a satisfy 

(5.35) 0<A o <2^-^. 

If€, K, K np satisfy 

2 — A a 1 2£ £ - 1 K 

(5 - 36) 2^3^ < * < d ' 2^V < ^ < V' K "*<-c 

then 



(5.37) 



AF < 
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Proof. The condition on K comes from (|5.9[) and Q5.27[l 9. The interval where 
K lies is not empty if A a < | and £ > 2 2 T 3 ^° ; together with (|5.27|b this gives the 
assumption required on £. In turn, it is possible to choose £ in such interval if (|5 .35[) 
holds. Remark that 2^5^ < § 5 so the previous condition on A a holds. Therefore the 
assumptions of Propositions [531 15.101 and !5.12l hold and then (|5.37[) follows. □ 

Proposition 5.14 (Global decreasing). Let to > 0; assume 115.35)) , H5.36)) , 

(5.38) L(Q+) < ^ 

and that the approximate solution U is defined in [0, T]. Then L(t) < to for any 
t G [0, T]; as a consequence, condition 15.19)) holds for any 1 or 3 wave in U . Finally, 
AF(t) < for all times t G [0, T] . 

Proof. Because of Proposition 1 5 . 1 31 for any time t not of interaction we have 

L{t) < L e (t) < F{t) < F(0+) < £L(0+) + KQ{0+) . 

Moreover Q(0+) < L(0+)L cd < L(0+)A o ; thanks to (|Q8|) and ([Q6]) we have 

(5.39) L(t) < L(0+)-(C + KA o ) < L(0+)-(2^-l) < m. 

Therefore the size |e| of any wave of families 1, 3 at time t > satisfies |e| < to, so 
1(05)1 holds. □ 

Remark 5.15. From (|5.38[1 we see that, in order to have L(t) < to, the smaller 
is £ the larger can be chosen L(0+). 

6. The convergence and the consistence of the algorithm. In this section 
we prove Theorem 12.21 We show first that for fixed v the algorithm introduced in 
Section |4] gives an approximate solution defined for every t > 0; more precisely we 
prove that at every time the number of interactions is bounded. Then we prove 
that the total amount of non-physical waves in each approximate solution is very 
small. The convergence of a suitable subsequence is assured by Helly's theorem; then 
consistence follows. 

6.1. Control of the number of interactions. We prove first that the size of 
the rarefactions in the scheme is small. 

Lemma 6.1. Consider a rarefaction of size e; then 

(6.1) |£|<7?e^. 



Proof. We analyze all possible situations. When the rarefaction is generated, one 
has < e < r\. When it interacts with a 1- or 3-wave, the size does not increase, 
see Remark 15.81 By Proposition 15 . 121 the size does not change when interactions with 
non-physical waves occur. 

The last case to be considered is when a rarefaction interacts with a 2-wave. In 
this case the size may increase; however, a rarefaction can meet a fixed 2-wave only 
once. Consider the case of a 1-rarefaction of size 6%, as in Proposition 15. 3i the other 
being analogous. If S 2 < then the size decreases, see (|5.13|) . If 62 > by (|5.12|) we 
have 

\ei\ = \6 1 \ + \e 3 \<\5 1 \(l + ±\6 2 \ \ < . 
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Summarizing the three cases above, we get |e| < ?]e L i<i/ 2 (or |e| < 7ye L <=<i/ 2 ) for a 1- 
rarefaction (resp. 3-rarefaction), where L^ d is the sum of the 2- waves at the right or 
left of the rarefaction. Then (|6.ip follows. □ 

Next, we prove that the number of interactions remains bounded in finite time, 
so that the approximate solution is well defined for all t > 0. We give first a lemma. 

Lemma 6.2. Consider the wave-front tracking algorithm described in Section^ 
under the assumptions of Proposition \5.14\ Then 

(i) the number of interactions involving a 2-wave and solved by the accurate 
Riemann solver is finite; 

(ii) the number of interactions where a new rarefaction of size s > rj arises is 

finite. 

Proof. Consider first (i) and refer to Proposition [531 Then, using (|5.16j) . we have 
< p(£ + K(A a — 2)/2) < 0, hence F decreases by a uniform positive quantity; 
since it is non increasing, this can happen only a finite number of times. 

Consider then (ii). After (i), it remains only to consider the case of two shocks 
of the same family interacting. 

Under the notation of the corresponding case in the proof of Proposition 15 . 101 we 
have e = e\ > r\ and 

AF < |ei| (1 - f + KA ) < ry (1 - £ + KA ) < 

because of (|5.36j) . Arguing as in (i), this can happen only a finite number of times. □ 
About (ii) in Lemma l6.2[ recall that if e > rj then the new rarefaction must be 

split into more than one wave. Therefore Lemma 16.21 can be rephrased by saying 

that, except for finite interactions, the number of waves emitted in an interaction is 

at most three, and this case occurs precisely when a non-physical wave is generated; 

moreover, in every interaction at most one wave per family is emitted. 

In a schematic way, apart from a finite number of interactions, in our algorithm 

the following holds (we will consider the set of non-physical waves as a 4 th family of 

waves) : 

(a) the interaction of an i-wave, i = 1,3, with a 2-wave is solved by a single 
i-wave, a 2-wave and a 4- wave; 

(b) in the interaction of just 1- and/or 3-waves, there is at most one outgoing 
wave of each family 1 and 3; 

(c) the interaction of a i-wave, i — 1,2,3, with a 4- wave is solved by an i-wave 
and a 4- wave. 

The following proposition is inspired by [3l Lemma 2.5]. 

Proposition 6.3. Consider the wave-front algorithm described in Section^and 
assume in the strip [0, T) x R the following: 

for some a\ < a2 < < b\ < 62 the waves of the first (third) family 

have speeds in the interval [01,02] (resp. [bt, 62]^- 
Then the number of interactions in the region [0, T) x R is finite. 

Proof. Assume by contradiction that in the region [0,T) x R there exists an 
infinite number of interactions. Unless of taking a smaller T we can assume that the 
number of interactions is finite in every strip [0, t] x R, < t < T, and that T is an 
accumulation point for the times of interaction. Then there exists a sequence (tj, Xj), 
j = 1,2,... of interaction points such that 

< tj < tj + i < T for all j and {tj,xj) — > (T, x) 
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for some x. Denote J — {(tj,xj):j = 1,2,...} the set of all interaction points; the 
set J is bounded in the strip [0, T) x M because of the finite propagation speed. 

The situation described in items (a)-(c) above holds except a finite number of 
interactions; let r < T be the maximum time of these "exceptional" interactions. It 

is not restrictive to assume that r <tj < T for all j = 1, 2 

Starting from a point of J~, we "trace back" all the segments up to t = 0; we 
repeat the procedure for all the points of J and call J- the set of the "traced segments" 
obtained in this way. In other words, a segment belongs to the set T iff it can be 
joined forward in time to some point of J by a continuous path along the wave fronts. 
The set T is not empty: for instance, two segments interacting at the point (tj,Xj) 
belong to T, for any j = 1,2, . . .. Observe the following dichotomy property of T 
which is used just below: 

two interacting waves either belong both to T or none of them does; 

moreover if at least one of the outgoing waves belong to J 7 , then both 

the incoming waves must belong to T . 
We partition now all the interaction points of the algorithm that occur for times t > t 
into the following sets: 

• Xq: the interaction points where no ingoing wave belongs to T\ 

• X\: the interaction points where both incoming waves belong to T and at 
most one outgoing segment belongs to T\ 

• Z2 : the interaction points where exactly two outgoing segments belong to T; 

• X 3 : the interaction points where three outgoing waves all belong to T . 

Because of the dichotomy property quoted above no outgoing wave in case Tq can 
belong to T . On the contrary, both incoming waves in Xi, X%, X3 must belong to T . 

Recall that we are considering times t > r. Therefore the maximum number 
of emitted waves in an interaction is three, and this happens only in the situation 
considered in I3, that is for interactions as in (a) above. The case of more than one 
emitted wave per family cannot occur, and so the outgoing waves in X2 belong to 
different families. By definition we have J n 2"o = and so J C X\ U Z2 U I3. 

Let V(t) be the total number of wave- fronts of the families 1, 2 and 3 that belong 
to T, at time t. The functional V(t) is non-increasing, and it decreases at least by 1 
across X\. Then X\ is finite. As a consequence, all the interaction points of J belong 
to I2 UI3, except at most a finite number. Let t\ 6 [t,T) be a time such that all 
points in X\ lies in t < t\. 

Let P — {x\, . . . jXjvj} the set of points of the x-axis where a 2-wave is located, 
Nx < N D . We consider two cases. 



x £ P . In this case we can choose the time t\ < T such that after that time 
no segment belonging to T crosses a 2-wave. Then all the points in J with tj > t\ 
belong to Xi- The same argument of [31 Lemma 2.5] can now be used, reaching a 
contradiction. 
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Consider (t*,x*) S J with t* > t\. As in [3J Lemma 2.5] we define for 
t G [t* , T) two continuous paths "fe(t), J r (t) starting at (t*,x*) in the following way. 

At (t*,x*) there are cither two or three outgoing segments belonging to T\ for 
times t > t* and sufficiently close to t* wc define 7^(i) to be the segment on the left 
and 7r(£) the one on the right. When ji(t) ("f r (t)) arrives at an interaction point, 
it is prolonged by a segment of T\ since they are at least two, it follows the one on 
the left (resp., on the right). Then the path 7^(i) is made by segment of the families 
1, 2 or 3, while j r (t) is made by segments of families 3 or 4; in fact a look to the 
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proof of Proposition 15 . 1 21 shows that the interaction of a 1-wave with a 2-wave always 
produces a 4-wave. 

We claim that the speeds of the paths 7^(f), r y r (t) are strictly separated. This is 
clear for a finite number of nodes among segments. If it happens that 7 r (t) follows 
once a 4-segment, then from that point on it will follow only 4-segments, so 72(f) < 
62 < s — j r (t)- If on the other hand "f r (t) never follows a 4-segment, then 7^(t) < < 

h < 7r(t). 

This proves the claim; the same argument of [3] then leads to a contradiction. □ 

6.2. Control of the total size of non-physical fronts. Assume as above 
that the assumptions of Proposition 15.141 hold. We assign inductively to each wave 
a a generation order k a as in [21 p. 140]. This is done according to the following 
procedure. First, at time t = each wave has order 1. Second, assume that two waves 
a and j3 interact at time i; if a and j3 belong to different families, the outgoing waves 
of those families keep the order of the incoming waves, the other waves assume order 
max{fc Q , kp}+l; if a and (3 belong to the same family, the outgoing wave of that family 
takes the order min{fc Q , fc^}, the other waves are assigned order max{fc Q , kp\ + 1. 

When specialized to the current setting this has the following consequences: 

• every 2-wave has order 1; when a i-wave, i = 1,3, of order fc interacts with 
a 2-wave, the outgoing i-wave has order fc, the other outgoing wave (of the 
family j, j = 1, 3, j ^ i, or a non-physical wave) has order fc + 1; 

• in the interaction of a 1- with a 3-wave the waves cross without changing 
order; in the interaction of two waves a, (3 of the same family i = 1,3, the 
outgoing wave of the family i takes order min{fc a , kp}, the wave of the family 
j = 1, 3, j =/= i, has order max{fc a , hp} + 1; 

• when a non-physical wave interacts with any other wave, both waves cross 
without changing order; in particular a non-physical wave keeps the order it 
has been assigned when generated. 

For t > not an interaction time and any k = 1, 2, . . . define, see (|5.2I) . (15.11) . 

V k (t)= ]T | 7 |+££ \l\+K np £ |7l 

•y>0 " y< °, "yGA/"-P 

Q k (t)= J2 E IfclM 

73 at the left of 5 2 71 at the right of 5 2 

F k (t) = V k (t)+KQ k (t) 

and 

^k(t) = X)VKt), Qk(*)=EQ/(f)> F k (t) = V k (t)+KQ k (t). 

t>k £>k 

We remark that 

F 1 (0+) = L e (0+) + KQ(0+), F k (0+) = for fc > 2. 

Observe that if a non-physical front interacts with another wave, the functionals above 
do not change; the same holds for interactions between 3- and 1-waves. Then we focus 
on interactions of waves of the same i family, i = 1, 3 (as usual denote j = 1, 3, j ^ i) 
and on interactions between 1- or 3-waves with a 2-wave. 

For h € N, denote by Ih the set of times t when an interaction occurs between 
two waves a and (3 of families 1 or 3 with maxjfccfc^} = h; denote by Jh the set 
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of interaction times t of a 1- or 3-wave of order h with a 2-wave. Finally, denote 
T h = I h U J h and / = U ft >i 4, J = U>i Jh, T = IU J. 

In order to control the total size of non-physical fronts we must strengthen the 
assumptions (|5.35[) and (|5.36p required in Proposition ^. 14l First, for any fixed m > 0, 
instead of (|5.35|) we require the stronger condition 

(6.2) < ,1„. ' 



2- Vd 



Then denote 



(6.3) A = , A 2 = -— -r-r -, ^=max^A,A 2 , 



£ ' K(2-A a )-C n K 

We need < /i < 1. From (|5.36[) we have K np C /K < 1; moreover we have < A < 1 
and A 2 > because of (|5.36jb . At last A2 < 1 holds iff 1 J A < K; this condition is 
stronger than the left inequality in (15.36|) t. Hence, instead of (I5.36|l we assume 



Remark the new upper bound required on £. As in the proof of Proposition 15.151 the 
interval where K varies is not empty if A Q < i and £ > ^Z%a • ^ n turn, we can find 

£ satisfying if (pT2"|) holds; remark that < ±. The last condition in (|6~4")l 

coincides with that in (|5.36|) . 

Remark 6.4. Proposition 15. 141 still holds under the stronger assumptions (|6.2[) . 
(16. 41) under the same condition (|5.38j) . because the inequality on the right hand side 
in (|6.4[) 9 has not changed, see (|5.39[) . 

Proposition 6.5. Fix m > and assume I6'.g|) , |6'.^| ). We Ziawe: 

1. r &T h ,h<k~ 2: then AF k = AF k = . 

2. t G T k -i: then AF k -i < 0, AF k = AF k > and 



fc-2 



(6.5) [AF k ] < (i([AF k -i]_ - [A*i] 



3. r G 7h, h > k: if h — k then AF k < 0; m any case AF k < and 

fc-i 



(6.6) £ l AF *L < t A ^L ■ 



£=1 



Proof. As we pointed out above, for r S 7 only interactions of waves of the same 
family are taken into account. Remark now that for interactions between 3-waves 
we have AQ e = L+ d AV e while AQ e = L~ d AVi holds for 1 -waves. So if t G I then 
AV t > (AVe < O^AVg = 0) implies AqI> (resp. AQ e < 0, AQ e = 0). 

Remark also that by Proposition l5.14l 

fc-i 

(6.7) J2 AF ? + Ap k<0- 
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OLh Pi 

Fig. 6.1: Interactions of 3-waves; h > I denote generation orders. 

1. If h < k — 2, no waves with order > fc are involved and then AFk = AFk = . 

2. Let h — k — 1. First, consider t g Ik-i', then AVk = AVk > 0. We prove 

that 

fc-2 

(6.8) [AV>,] + < -([AV*_ X ]_ - £ t A ^]+) • 
Indeed, from ()5.30|) - (|5.32p . we deduce 

fc-2 

(6.9) e[AVfe] + + AVfc_i+^A^<0. 

If min{fc Q , A;,g} = fc-1, then AVi = for t = 1, . . ., fc-2 and ^6T9j becomes £[AVfc] + + 
AVfc-i < 0; this implies AV fe _i < and hence [AVfc]+ < (l/^)[A^_i]_, that is ([H78]> . 
If minjfcQ,, fc^} = £ < fc — 2, then AVfc_i < since no waves of order k — 1 are present 
after the interaction; therefore the estimate (16. 9j) becomes 

£[AF fc ] + - [AV fc _i]_ + AV/ < . 

If AV e > we get (fTS}. If AV^ < we check directly that £[AVfc]+ - [AV fe _i]_ < 
^ 2 |<sri I — I o^3 1 < because of (|6.4lh . This proves f|6 .8[) : from this proof we see moreover 
that AVk-i < and then AF fe _i < 0. 
From ([678]) and < AQ k < A AV k we have 

/ fc ~ 2 \ 

(6.10) < AF k < (1 + if A )[AT4] + < A( [AT4_i]_ - £ [Ay £ ] + 1 . 

^ e=i ' 

We prove now that 

fc-2 

(6.11) [AQ fc _i]_ - ^ [AQ £ ] + >0. 

£=1 

We have only to consider the case in which [AQ£] + > for some I < k — 2; but in 

this case [AQ fe _i]_ - EiZii^Qeh = ^(M + |#a| - N) > because of (jOU). 
Therefore ([631) for r € J fe _i follows from (I6T0]) and l|H3T]) . 

Second, assume t 6 Jfe_i; we prove that 

(6.12) [AF fe ] + < /i[AF fe _i]_ . 

22 



Indeed, if the reflected wave is a physical wave then, under the notation of Proposi- 
tion [Ol 

< g , A<y fe _i < — |oxd 2 1 H — A = — (2 - A ) 

so that AiVi < [f - i<T(2 - A )] \5i5 2 \/2 < because of ([6~4]> . Then (j6T2j) follows 
since = AF& > and 

[AF,] + < Ml(£ +i ^ o ) - [ X (2 - A>) - £]-A a < A 2 [Afn]_ < ^[A^_ 1 ]_ . 

If the reflected wave is a non-physical wave we have, under the notation of Proposi- 
tion EH 

< AF k = AV k < K np C \S6 2 \ , AV fc _i = , AQ fe _! = -|<« 2 | 
and then 

[Afn] = #|<tt 2 |, [AF k ] <%^[AFn] < M [AFn] . 

The estimate (|6. 12[) is then completely proved. From (|6.12|) we get (|6.5|) since no 
waves of order < k — 2 are involved in the interaction. 

3. Finally, let h > k. We first consider t £ 1^, h > k. If min{fc Q , kp} > k, 
then AVfe = AL^ < 0. If min{fc a , kp} < k — 1, assume k a > k and kp < k — I; then 
A14 < £M - M < l)|a 3 | < by (jOSjh . This proves AVJ. < and then 

AF fc <0. 

From (|6.7|) we deduce that J^=i A-Fi < [A-Ffc]_. Since at most one non-zero term is 
present in the first sum, (|6.6p follows. In the case h — k we have A 14 = [AVfc+i]+ + 
AV k < 0; hence AV k < and then AF k < 0. 

Now assume t £ Jh, h > k. Then AF k — AF < and no waves of order < k are 
present, so (I6.6P holds. If /i = k and the reflected wave is physical, from the proof of 
Proposition 15.31 and (|6.4[) we find that 

AF k = AT/ fe + XAQ fc < (£ - + KA ) < . 

If h = k and the reflected wave is non-physical then AV k = 0, AQ k < and AF k < 0. 
□ 

Summarizing, for r e T we have the following table: 



%\ h 


<k-2 k-1 k >k+l 


AF k 


+ - ± 


AF k 


+ 



We write Fjr(t) = ^2 T < t [AFk(r)]± for k>2. For simplicity the time r in such sums 
is omitted. 

Lemma 6.6. Under the assumptions of Proposition HT5l we have 

(6.13) F+(t) < y. (L € (0) + iTQ(0)) + ]T [ Ai "i] + 

fc-2 

(6.14) F+(t) < x (t) + £ [AF fc _ x ] + -EE t A ^]+) ' ^ 3 ■ 
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Proof. By Proposition 16.51 see also the table above, the functional increases 
at times r € Tk-i, decreases at times r € 7^, while it does not have a given sign at 
times t <E Th, with h > k + 1. 

First, by summing up (16. 5|) we obtain 



fe-2 

(6.15) F+(t) < M £ ([AF t -i]_ - J] [AJk]. 

for fc > 2, where the last term in (|6.15p is missing if k = 2. 
Recall now that Fi(0) = L c (0) + ifQ(O); therefore 

Fi(t)<L i (0) + KQ(0)-Y,[AFi]_+ ]T [AF X ]_ 

and then 

(6.16) J2^L< L dO)+KQ(0)+ ]T t AF i]+- 



On the other hand Ffe(O) = for k > 2; from Proposition 16.51 we have 
F fe (i)<£[AF fc ] + -]T[AF fc ]_ + £ [AF fc ] + . 

T k - ± T k T h ,h>k+1 



Moreover 



and then 



y £[AF k ] + = y £[AF k ] + = Fjt(t) 



(6.17) £[AF fc ]_<F+(i) + 2 [ AFfc ]+- 

T fc T h ,h>k+1 

From (j6~T5l) . (|636jl . ([6~T7f we get (j6~T3| . ([614]) . □ 

Proposition 6.7 (A contraction property). Under the assumptions of Proposi- 
tion \6.5[ for any t > and fc > 1 we ftcwe 

(6.18) V fc (t) < F k (t) < /Z" 1 ■ (L 5 (0) + A-Q(O)) . 



Proo/. The estimate (|6TT5)| holds for fc = 1 because = F(t) < L^(0)+KQ(0). 
Next we prove by induction on fc > 2 that for any t 

fe-i 

(6.19) F+(t)<^(L i (0) + KQ(0))+ £ ^ [AF £ ] + . 

T h ,h>k 1=1 



Since by summing up (|6 . 6[) we obtain 

fc— l 

(6.20) F fc -(t)> £ Et A ^ + 

T h , ft>fc £=1 
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then (|5TT5j) will follow from (j6~E))) for any fc > 2 because of 

Formula (j6~IT>)) for fc = 2 reduces to (|6"7T5|) . Next, assume that (|6.19p holds for 
some fc > 2. By (|6.14p and the induction assumption 

r h ,h>fc+i r fc £=i 

fe-i fe-i 

< M fc (L C (o)+xQ(o)) + / i( £ Et A ^] + + E i A ^-] + -EE[ A ^] + ) 

T h ,h>k£=l T h ,h>k+1 T fc «=1 

fc 

</(^(0)+ifQ(0)) + M ^ Et A ^] + - 

T h ,h>k+1 1=1 

Since /x < 1, we get (pl^) for fc + 1. □ 

Remark 6.8. We comment now the case A a = 0. In this case system (jl.lj) 
reduces to the p-system with pressure law given by (12. for fixed A. According to 
our front-tracking algorithm, stationary and non-physical waves do not appear and so 
Led — Q = 0; the algorithm reduces to the one introduced in [3]. Then, Proposition 
IG. 51 holds with Vk and l/£ replacing Fk and /u, respectively, and at last (|6 . 1 8[) reads 

V k (t)< _L.L f (Q). 



Next, we conclude that the total strength of all non-physical waves is small by 
proceeding as in [9l page 142]. Recall the notation in Section 01 First, using (|5.39| . 
we find that the sequence {v y } related to the initial data (v^,u^,X^) is uniformly 
bounded and in particular uniformly bounded away from 0; then the eigenvalues e\ 
and e3 are bounded and this makes possible the choice of a suitable s. We have two 
more parameters 77, p to be chosen. Fix 77 > with the condition 77 = r\ v — > as 
v — > 00 and estimate the total number of waves of order < fc. We have 

E M(*)<v*(t) + ]T | 7 |(*) 

7GJVP yeNP,k y <k 

<^-(^(0)+if(5(0)) + C p ■ [number of fronts of order < fc] < - 

by choosing fc sufficiently large to have the first term < 1/(2^) and then choose p 
small enough to have the second term < 1/ (2v) . 



We accomplish now the proof of Theorem 12.21 Define first 



(6.21) fc M = i^p. 

From the properties of the function d(m) stated in Remark l5.7l we see that fc(0) = 1/2 
and that kirn) is decreasing, tending to for m — > +00. The assumption (|2.8p implies 
that (|6~2| holds. 

Now, by hypotheses (|2.7|) it follows that we can choose £ such that 



I 1 771 

2 TVl0g ^) + 2^ TV ^) < —1 < {l ~ 2A ° )m 
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and that (|6.4p i holds. Hence, using (|3.12p and (i) in Section [?J we have 

I 1 771 

L(0+) < -TV \og( Po ) + ——TV(u ) < 



2 °^ UJ 2infa v ' 2£ - 1 

so that the hypotheses (|5.38[) of Proposition 15.141 holds. Theorem 12.21 now follows 
along the lines of [9l §7.4]. 

Appendix A. The weighted total variation. In this Appendix we prove 
Proposition ^. II Remark that the map d(a, b) = ^7^ is a distance on R+, as one can 
easily prove. 

We start with the proof of the inequality on the right in (|2.5D . It is enough to 
prove that 

We claim that 

(A.2) logf > 2 ^~ 1 ' forOl 

K ' & ~ t+1 

where the inequality is strict if t > 1. To prove the claim it is sufficient to notice that 
the function 4>(t) = logi — 2|^j- vanishes in 1 and 4>'(t) — jtj^fp > if t > 1. 
We apply (|A.2[) to t = x/y for < y < x and arguing by symmetry deduce 

2\x — y\ 

I log a; — logyl > , for every x, y > . 

x + y 

Then (|A.1|) follows. The proof of the inequality on the left in (|2.5[) is analogous, 
starting from the inequality 

(A.3) -logt < 2(t ~ ^ forOl 

with strict inequality if t > 1. 

Assume now that / S C(R). To show that WTV(/) = TV (log(/)), we have to 
prove that the inequality 



(A.4) TV(log(/)) <2sup]T 



\f( Xj ) - /(^-_x)| 



holds for any interval [a, 6] C K. Consider any partition {x Dl x%, . . . , x n } of the interval 
[a, 6]. By the mean value theorem for some Q between f(xj) and f(xj-i), and by the 
intermediate value theorem applied to /, we get 

|log (/(*;))- log (/fe-r))| = 

Vim) + 

for some r\j € [xj-i, Xj]. We exploit again the continuity of / in [a, b]. On one hand its 
image is compact, then min[ a u f = m > 0. On the other / is uniformly continuous 
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in [a, b], so that for any e > there exists S e > such that |/(a;) — f(y)\ < e if 
\x - y\ < S e , for x,y S [a, 6]. 

Fix now any e > 0; without loss of generality we can consider partitions of the 
interval [a, b] of mesh less than 5 e . Assume for instance f(xj^i) < f(xj); then from 
the inequalities 

f( Xj )-f( Xj -i) < e, f(xj-i) < f(vj) < f(xj) 

it follows 

f(Xj) + fjXj-l) < 2/(Xj_i) +£ £ 

2/(^) " 2/(x i _i) " 

The inequality (|A.4jl follows by remarking that then for any partition of mesh less 
than S e 

± Hog u M) - io g (/fe-,,)! < (i + i) 2 g ^fgf ■ 

The proof of Proposition ^. H is complete. 

Remark A.l. Observe that, if TV(log(/)) < oo and / is discontinuous, then 
the inequality on the right in (|2.5[) is strict, because of the strict inequality in (|A.2() . 
For example, if / has a single jump and assumes the values c > and d > 0, then 
WTV(J) = < | logc - logrf| = TV(log(/)). 

Remark moreover that WTV and TV are not equivalent, in the sense that there 
does not exist a positive constant C such that C ■ TV(log(/)) < WTV(/). This 
follows from the fact that clearly the inequality Clogt < does not hold for 

every t > 1. 

Appendix B. Shock-rarefaction interactions. In this Appendix we consider 
a particular case of Lemma 15.61 Actually, this is the only case needed in order to 
define a decreasing functional, Section [51 however, we needed further analysis for the 
control and treatment of the non-physical waves. 

Lemma B.l (The case SR, RS — * SS). Consider the interaction of a shock &i 
and a rarefaction Pi of the same family, i = 1,3, producing two outgoing shocks S\, £3. 
Then there exists a smooth function B satisfying |otj| < B(cti) < min{sinh(|a;|), 2|aj|} 
such that 

(B.l) < fa < B{oti) . 

Moreover, assume 

(B.2) < m 

for some m > and denote c = c(m) = cos |^ m ) 1 Then both the variation of shock 

J V / cosn(m) + l J 

waves and the reflected wave Sj, j ^ i, j = 1,3, are estimated by the interacting 
rarefaction as 

(B.3) M + N-N < (2c- 1)- 1/3(1, 

(B.4) | ei |<c.|A|. 
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Proof. We focus on the case i = 3, j = 1, see Figure [?75f bV Therefore we consider 
a 3 < 0, fi 3 > and e x < 0, e 3 < 0. Then ([5T5J) . ([5TB]> become 

(B.5) |ei|-N= l&| -W 

(B.6) sinh(| ei |)+sinh(|e 3 |) = sinh(|a 3 |) - |/3 3 | . 

From the second equation |ex| < |«3 1 , | £3 ] < |a 3 | and |/3 3 | < sinh(|a 3 |); using the first 
equation and |e 3 | < |a 3 |, we get |ei| < |/3 3 |. Therefore in conclusion 

(B.7) | ei | <min{|a 3 |,|/3 3 |}, |e 3 | < |a 3 |, |/3 3 | < sinh(|a 3 |) . 

Step 1: notation. We set x = \/3 3 \, y = \ei\, z = |a 3 |, so that 
(B.8) |e 3 | =y-x + z. 

Under this notation, (|B.6p writes as 

(B.9) F(x,y; z) = sinh y + sinh(y — x + z) — sinh z + x = 

for x > 0, y > 0, z > 0, y — x + z > 0. By (IB.7j) . any solution of (IB.9|) satisfies 

(B.10) y<z, y<x, x < sinh(^) . 

y - x + z = |e 3 | 




Fig. B.l: Interactions. 

Step 2: the threshold. Observe that, despite the last inequality in (|B.7|) . we may 
well have |/3 3 | > \a 3 \, that is, x > z. Consider in fact the limit case of e 3 = 0: we 
have y — x — z > and sinh(y) = sinh(z) — x that give 

sinh(a; — z) = sinh(z) — x, x > z . 

The last equality is the relation needed for /3 3 , a 3 in order to have that the shock 
and rarefaction cancel out exactly, giving rise only to a wave of the opposite family. 
Observe that the size of the rarefaction must be larger than the one of the shock. 
Under the notations above, the threshold curve separating the case of the outgoing 
waves SiS 3 from the case S\R 3 is given by 

(B.ll) f(x, z) = sinh(x — z) — sinh(z) + x = . 

Since f x — cosh(x — z) + 1 > and /(z, z) < 0, the implicit equation f(x, z) = 
is solved by x — x a (z) > z with x' a {z) — ^^^l^yl^^ > for every z > 0; the 
curve has for tangent at (0,0) the line z = x. Observe that x (z) < 2z because 
/(2z, z) = z > 0. In conclusion 

(B.12) z<x (z)<2z. 
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z = \a 3 \ x = z 

o o -f ^ n / threshold x = x (z) 



y = x 




x = \fo\ 



(a) 




x (z)/2 z x o{z) x 
(b) 



Fig. B.2: (a), the threshold curve f(x,z) — sinh(a; — z) — sinh(z) + x = 0; 
(b), the domain D z and the function y — y(x; z). 



This estimate and the third inequality in (|B.7|) prove (|B.1|) for B{oti) = x (\ai 
can prove more than (|B.12|) . that is, 



We 



(B.13) 



lim (x (z) - 2z) = . 

2 — >+00 



Indeed we show that the inequality x (z) > 2z — q holds for large z and q > . This 
follows from f(2z — q, z) = sinh(z — q) — sinh z + 2z — q ~ e z (e -9 — 1) /2 — * — oo for 
z — * oo. 

Remark that the fact that £3 is a shock implies that 
(B.14) f(x, z) = sinh(a; — z) — sinh(z) + x < . 

Step 3: the amount of shocks can increase. From (|B.5j) we have 
(B.15) |ei|+N-|a3| = 2|ei|-|ft|. 

We prove now that the inequality |ei| + |e3 j — \a 3 \ < 0, or equivalently |ei| < i | /3 3 1 , 
does not hold if m is large. 

The equation giving |ei| = y in terms of — x, for a given parameter \as\ = z, 
is (|B.9|) . to be considered in the domain 



{(x, y): < x < x {z) , y > max{0, x — z}} 



for z > 



where x (z) satisfies the equation (jB.llj) . Since F y = coshy + cosh(y — x + z) > 0, 
the implicit equation (|B.9j) defines a function y = y(x; z) with < y(x; z) < x 
and y(x;z) < z, see (IB. 101) . Remark that F(x,x;z) = sinh.T + x > 0. Moreover 
F(x, 0; z) — sinh(z — x) — sinh z + x so that F^O, 0; z) — 0; by F x = 1 — cosh(z — x) < 
0, we deduce F(x,0;z) < if x > 0. By the implicit function theorem we have 
y'{x;z)>0, 



(B.16) 



, . , cosh(z) — 1 
l/(0;z)= , ,, €[0,1) 



cosh(z) + 1 



and y"(0;z) = --^ 



4 sinh(z) 



< 0. The function "^l^ll , \ is increasing and then 
for z S [0, m] its maximum is e°sh(mj+i ' ^ n ^ s quantity is strictly larger than 1/2 if 



sh(z)) a 



sh(z) + l 



cosh(m) — 1 
cosh(m) + l 

m > log(3 + 2y2)- Thus in general the estimate y(x; z) < x/2 cannot hold 
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Step 4: proof of the estimate. From (|B.15[) we see that |ex | + |es| — 1 0:3 j < 
(2c - 1) • \/3 3 \ |ei| < c • \0 3 \, that is, (|B.3|) and (|B.4[) are equivalent; we shall 

prove (|B.4p . To bypass the study of the function y{x\ z) we define 

<f>(x; z,c) — F (x, cx; z) = sinhfca;) + sinh (z — (1 — c)x) — sinhz + x . 

If 1/2 < c < 1 then z > (1 - c)x and $(0;z,c) = 0, 

$ x (%', z, c) = 1 + ccosh(cx) — (1 — c) cosh (2; — (1 — c)x) , 
$ xx (a;; z, c) = c 2 sinh(ca;) + (1 — c) 2 sinh (z — (1 — c)x) > . 

Therefore the function x — > <fr x (x; z, c) is increasing and then z, c) > if 

•^(O; z, c) = (1 + c) - (1 - c) cosh(z) > , 

that is if cosh(z) < this is just (|B.2[) . Then y(x;z) < cx for all a; S (0, x (z)) 
and so (|B.4p is proved. □ 

Remark B.2. From (|B.16j) we see that condition (|B.2j) is equivalent to the 
geometric condition y'(0;z) = (coshz — l)/(coshz + 1) < c. Moreover, as we noticed 
in the above proof, condition (|B.2[) is equivalent to 



(B.17) cosh(H) < i±£ 



which, in turn, is equivalent to 



(B.18) K-|<log 



1 



From the definition of the strength, one has that \ai\ = (l/2)log(v max /v m i n ), where 
v-max — max{wf,w r }, v m i n — min{t^,tv}, being vi, v r , respectively, the left and right 
values of v for the wave of size cti. Hence (|B.18|) is equivalent to 



'max ^ (1 + ^) 



1 - C 



Remark B.3. In the proof above we showed that AL shocks = |ei| + |£a| — 
I CK3 1 = 2|ei| — |^3| may be positive, differently from Nishida's paper, where it is 
always decreasing. This depends on the definition of the wave strengths, which was 
imposed to us in order to have good estimates when dealing with interactions with 
the 2 waves. In any case AL = AL shocks + AL rarcfac tions < 0. 

Remark B.4. Under the notation and assumptions of Lemma fB. II we verify that 

(B.19) \ej\Kc\oi\. 

This estimate, together with (|B.4j) . allows us to obtain (in a special case) the analogue 
of (|5.20p with c in place of d. 

The proof makes use of a numerical computation. As in that lemma we consider 
the case j = 1, i = 3. Let z > be fixed and < x < x (z). From the proof of 
Lemma iB.ll we deduce y — y(x, z); for z fixed the function x y{x, z) is increasing. 
Define then 



Y(z) = y(x (z),z) = x {z) - z . 
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In order to prove (|B.19p it is sufficient to prove that 



Y(z) < cz if cosh(z) < . 

1 — c 

Remark that from (|B.13|) we know that Y(z)jz — ► 1 for z — > +00; however the 
constraint (|B.I7[) implies that z is bounded. The inequality Y(z) < cz is equivalent 
to x {z) — z < cz, i.e., x (z) < (1 + c)z. Therefore we need to prove that 

1 + c 

(B.20) cf)(z, c) = sinh(cz) - sinh(z) + (1 + c)z > if cosh(z) < . 



Notice that cosh(z) < means < z < z c for z c = log(f±§ + J (y^f) 2 - 1). 

Formula (|B .20|) is shown to hold true by numerical computations. Remark however 
that 

(p'(z, c) = ccosh(cz) — cosh(z) + 1 + c 
4>"(z, c) = c 2 sinh(cz) - sinh(z) < 

so </>'(0,c) — 2c, (f>(-,c) is concave, lim z ^ +oc </>(z, c) = —00 and </>(•, c) has a single 
point of maximum; at last <j)'{z c ) — c (cosh(cz c ) — cosh(z c )) < 0. This concludes the 
proof of dHUSJ). 

Remark that if c = 1/2 then (f>(z, c) = sinh(z/2) — sinhz + |z. If z is such that 
coshz = 3, then sinhz = 2\/2, sinh(z/2) = 1, z = log (3 + V$>), and (|B.20|) holds with 
strict inequality. 

Acknowledgments. The authors kindly acknowledge support by Progetto 
Gnampa 2005 Analisi Asintotica Per Sistemi Iperbolici Non Lineari. The authors 
thank Graziano Guerra for stimulating remarks and Umberto Massari for hints on 
BV functions. 

REFERENCES 

[1] R. Abeyaratne and J. K. Knowles, Kinetic relations and the propagation of phase boundaries 
in solids, Arch. Rational Mech. Anal., 114 (1991), pp. 119-154. 

[2] D. Amadori and A. Corli, A hyperbolic model of multi-phase flow. Proceedings of the Con- 
ference HYP2006, in print. 

[3] D. Amadori and G. Guerra, Global BV solutions and relaxation limit for a system of con- 
servation laws, Proc. Roy. Soc. Edinburgh Sect. A, 131 (2001), pp. 1-26. 

[4] F. Asakura, Wave-front tracking for the equations of isentropic gas dynamics, Quart. Appl. 
Math., 63 (2005), pp. 20-33. 

[5] , Wave-front tracking for the equations of non-isentropic gas dynamics, To appear, 

(2005). 

[6] S. Benzoni-Gavage, Analyse numerique des modeles hydrodynamiques d'ecoulements diphasi- 

ques instationnaires dans les reseaux de production petroliere, Ph.D. Thesis, Ecole Normale 

Superieure de Lyon, Lyon, 1991. 
[7] F. Bereux, E. Bonnetier, and P. G. LeFloch, Gas dynamics system: two special cases, 

SIAM J. Math. Anal., 28 (1997), pp. 499-515. 
[8] S. Bianchini, The semigroup generated by a Temple class system with non- convex flux function, 

Differential Integral Equations, 13 (2000), pp. 1529-1550. 
[9] A. Bressan, Hyperbolic systems of conservation laws. The one- dimensional Cauchy problem, 

Oxford University Press, 2000. 

31 



[10] A. CORLI AND H. Fan, The Riemann problem for reversible reactive flows with metastability, 
SIAM J. Appl. Math., 65 (2004/05), pp. 426-457. 

[11] C. M. DAFERMOS, Hyperbolic conservation laws in continuum physics, Gmndlchrcn dcr Math- 
ematischen Wissenschaften, 325, Springer- Verlag, Berlin, second ed., 2005. 

[12] R. J. DiPerna, Existence in the large for quasilinear hyperbolic conservation laws, Arch. Ra- 
tional Mcch. Anal., 52 (1973), pp. 244-257. 

[13] H. Fan, On a model of the dynamics of liquid/vapor phase transitions, SIAM J. Appl. Math., 
60 (2000), pp. 1270-1301. 

[14] J. Glimm, Solutions in the large for nonlinear hyperbolic systems of equations, Comm. Pure 
Appl. Math., 18 (1965), pp. 697-715. 

[15] L. GOSSE, Existence of L°° entropy solutions for a reacting Euler system, Port. Math. (N.S.), 
58 (2001), pp. 473-484. 

[16] T.-P. Liu, Initial- boundary value problems for gas dynamics, Arch. Rational Mech. Anal., 64 
(1977), pp. 137-168. 

[17] , Solutions in the large for the equations of nonisentropic gas dynamics, Indiana Univ. 

Math. J., 26 (1977), pp. 147-177. 
[18] Y. Lu, Hyperbolic conservation laws and the compensated compactness method, Chapman &: 

Hall/CRC Monographs and Surveys in Pure and Applied Mathematics, 128, Chapman & 

Hall/CRC, Boca Raton, FL, 2003. 
[19] T. NlSHlDA, Global solution for an initial boundary value problem of a quasilinear hyperbolic 

system, Proc. Japan Acad., 44 (1968), pp. 642-646. 
[20] T. NlSHlDA AND J. A. Smoller, Solutions in the large for some nonlinear hyperbolic conser- 
vation laws, Comm. Pure Appl. Math., 26 (1973), pp. 183-200. 
[21] Y.-J. Peng, Solutions faibles globales pour I'equation d'Euler d'un fluide compressible avec de 

grandes donnees initiales, Comm. Partial Differential Equations, 17 (1992), pp. 161-187. 
[22] , Solutions faibles globales pour un modele d'ecoulements diphasiques, Ann. Scuola Norm. 

Sup. Pisa CI. Sci. (4), 21 (1994), pp. 523-540. 
[23] F. Poupaud, M. RASCLE, and J. -P. Vila, Global solutions to the isothermal Euler-Poisson 

system with arbitrarily large data, J. Differential Equations, 123 (1995), pp. 93-121. 
[24] D. Serre, Systems of conservation laws. Vol. 1 and 2, Cambridge University Press, Cambridge, 

1999-2000. 

[25] J. Smoller, Shock waves and reaction- diffusion equations, Springer- Verlag, New York, sec- 
ond cd., 1994. 

[26] J. B. Temple, Solutions in the large for the nonlinear hyperbolic conservation laws of gas 

dynamics, J. Differential Equations, 41 (1981), pp. 96-161. 
[27] D. H. Wagner, Equivalence of the Euler and Lagrangian equations of gas dynamics for weak 

solutions, J. Differential Equations, 68 (1987), pp. 118-136. 



32 



